o 

(N 



SPIDER - V. Measuring Systematic Effects in Early- Type Galaxy Stellar 
Masses from Photometric SED Fitting 

R. Swindle^ and R. R. Gal 

Institute for Astronomy, University of Hawaii i 

F. La Barbera 

INAF - Osservatorio Astronomico di Capodimonte, Napoli, Italy 

R.R. de Carvalho 

CN ' Instituto Nacional de Pesquisas Espaciais/MCT, S. J. dos Campos, Brazil 

O ■ ABSTRACT 

u 

,J^ ] We present robust statistical estimates of the accuracy of early-type galaxy stellar 

masses derived from spectral energy distribution (SED) fitting as functions of various 
O I empirical and theoretical assumptions. Using large samples consisting of 40,000 galaxies 

"^ ' from the Sloan Digital Sky Survey, of which 5,000 are also in the UKIRT Infrared Deep 

^, ' Sky Survey, with spectroscopic redshifts in the range 0.05 < z < 0.095, we test the 

reliability of some commonly used stellar population models and extinction laws for 
computing stellar masses. Spectroscopic ages (t), metallicities (Z), and extinctions (A) 
are also computed from fits to SDSS spectra using various population models. These 
constraints are used in additional tests to estimate the systematic errors in the stellar 






•/^ I masses derived from SED fitting, where t, Z, and A are typically left as free parameters. 

C""^ ' We find reasonable agreement in mass estimates among stellar population models, with 

variation of the IMF and extinction law yielding systematic biases on the mass of nearly 
a factor of 2, in agreement with other studies. Removing the near-infrared bands changes 
the statistical bias in mass by only 0.06 dex, adding uncertainties of 0.1 dex at the 95% 

^ , CL. In contrast, we find that removing an ultraviolet band is more critical, introducing 

H I 2? uncertainties of 0.15 dex. Finally, we find that stellar masses are less affected by 

absence of metallicity and/or dust extinction knowledge. However, there is a definite 
systematic offset in the mass estimate when the stellar population age is unknown, up to 
a factor of 2.5 for very old (12 Gyr) stellar populations. We present the stellar masses 
for our sample, corrected for the measured systematic biases due to photometrically 
determined ages, finding that age errors produce lower stellar masses by 0.15 dex, with 
errors of 0.02 dex at the 95% CL for the median stellar age subsample. 

Subject headings: galaxies: elliptical, ETC, stellar mass - methods: SED fitting 
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Introduction 



Our understanding of the formation and evolution of early-type galaxies (ETGs) represents 
a key ingredient in models of galaxy formation. Here, we consider ETGs to be bulge dominated 
galaxies with passive spectra in their central regions. Observationally, they are characterized by 
elliptical isophotes, generally redder colors, and a sharp 4000A break, corresponding to an accumu- 
lation of absorption lines of mainly ionized metals, reflected in their rest frame UV/optical colors. 
This feature is typical of old stellar populations with little to no ongoing star formation. 

At higher redshift, most ETGs can only be studied in integrated light, and interpretation of 
their photometric and spectroscopic properties requires population synthesis models. These single 
stellar population (SSP) models usually consist of stars born at the sam e time with equal initial 



element compositions, evolved using the isochrone synthesis technique (see lCharlot fc Bruzuall 



1991 



for a review), where stars of different masses follow different evolutionary tracks. SSP models can be 
combined to produce arbitrary star-formation histories, although most studies restrict themselves 
to histories with exponentially declining star-formation. Such tau-models are parametrized by the 
e-folding time of this decline, r. The spectral energy distributions (SEDs) from a set of models 
with various parameters (e.g. initial mass function (IMF), star- formation history, age, metallicity, 
and extinction) is compared to the photometric or spectroscopic observations to derive a best-fit 
template. A fundamental parameter derived from such SED fitting is the overall normalization 
of the model relative to the observations, which gives the galaxy stellar mass content. Indeed, a 
measurement of the stellar mass of a n ETG is involved in many useful sca ling relations, such as the 
size-mass relation (|Shen et al.l l2003l ) and downsizing (jCowie et al.lll996l ). where the evolutionary 
history of ETGs is seen to follow different time scales as a function of their stellar mass content. 
Stellar mass assembly in galaxies is also used in tests of hierarchical models, such as the e volution of 



the number density and size of both early an d late-type galaxies as a function of redshift (JBell et al. 



2003 



Bundy et al. 



2005 



Ilbert et al.ll2010l ). Furthermore, upcoming surveys (e.g. PanSTARRS, 



LSST, DES) will provide only photometry, so it is crucial to understand how to obtain reliable 
stellar masses from SED fitting techniques. 

The degeneracies among the multiple model parameters which are required to reproduce the 
observed SEDs of the galaxies, along with the available photometric bandpasses, determine the 
uncertainty of the mass estimates. For example, models and obser vations show that the rest; 



frame near-infrare d (NIR) galaxy flux correlates well with stellar mass (JKauffmann Sz Chariot 



1998 



Cowie et al.lll996l ) due to weak contributions from hot, young stars and dust extinction at these 
wavelengths. But both different model predictions and/or absence of NIR data in fitting the 
SED ca n result in stellar masse s which differ by a factor of 2, for both low and high redshift 
samples (jvan der Wei et al.ll2006l ) . Of particular interest in this study, the effects of unknown stellar 
population age, metallicity, and extinction on the stellar masses derived for elliptic al galaxies have 
not been quantified. It is already well known that the age-metallicity degeneracy (| WorthevI 1 1 994 ) 
cannot easily be broken with broadband colors alone. 



All of these factors either depend on or affect the resu lting evolution of the spec t ral energy 
distrib ution , models of which have been gen erated by, e.g. iBruzual &: CharlotI (J2003l ). iMaraston 
(|2005l ) , and iFioc &: Rocca-Volmerangd (|l997l ) . Many phases in stellar evolution are still not well 
understood, but one key ingredient in these models, the thermally-pulsating asymptotic giant 
branch (TP-AGB) phase, is particularly influential for determining the galaxy stellar mass for a 
range of dominant stellar population ages. Light from these stars largely influence s the integrated 



bright ness of the NIR con tinuum, the effects of which were recently highlighted by IMaraston et al. 



(|200g) and iBruzuall (|2007j), who found systematic differences of factors of ~2 in mass and age over 
a large range in redshift. This example has shown us the effects of one sensitive parameter in the 
SED fitting process and demonstrates the need to further test the dependence of the stellar mass 
estimate on model parameters and the use of common ground-based observation filters. 



Several authors (e.g. IConroy et al.ll2009l : iTreu et al.ll2010l : iRaichoor et al.ll201ll ) have studied 
some of these effects on the photometrically derived stellar mass, quantifying systematic offsets 
between masses calculated against different model parameters. However, they all suffer from a lack 
of spectroscopic data which can provide independent constraints on many of the galaxy properties. 
Longhetti &: Saraccd ((20091) used a synthetic catalog of intermediate redshift {1 < z <2) ETGs to 
test the dependence of stellar population models and assumed model parameters on stellar mass. 
However, they compare a sy nthetic catalog to o nly 125 galaxies obs e rved i n the GOODS field, using 
a different SED fitting code. Illbert et al.l (|2010l ) and lBernardi et al.l (J2010l ) test the accuracy of large 
surveys (~10,000 and ~200,000 galaxies, respectively) at estimating photometric stellar masses 
compared to those from available spectroscopy. However, these projects select mostly galaxies at 
high redshift and neglect the possible combined effects of assuming various free parameters. 

This paper is part of a series examining the global and internal properties of ETGs in the 
nearby Universe, combining optical and NIR photometry with spectroscopic data. The Spheroids 
Panchromatic Inv e stigati on in Different Environmental Regions (SPIDER) project is described in 
La Barbera et al.l (|2010al . hereafter Paper I). The second and th ird papers of the series pr esent a 
thorough an alysis of the optical+NIR scaling relations of ETGs ( La Barbera et al.ll2010bl Ja). while 
in Paper IV (JLa Barbera et al.ll2010dl ) we have analyzed the optical+NIR internal color gradients of 
ETGs. Here we present an extensive comparison of stellar mass estimates obtained for low redshift 
ETGs by varying several observational and model parameters (presented individually in Section [7]) . 
The paper is organized as follows. Section [2] describes the galaxies, comprised of a carefully selected 
sample of high S/N, low redshift ETGs. In Section [3] we describe the different stellar population 
models used to compute model galaxy SEDs. In Section [5] we provide an overview of the general 
technique applied to fit the theoretical SEDs to the observed colors and the SDSS galaxy spectra. 
The reliability of these fits is assessed in Section [5l by comparing photometric and spectroscopic 
fitting results. Our method for determining stellar masses is described in Section [6l while the results 
of the various model comparisons are discussed in Section [Tl Throughout the paper, we adopt a 
cosmology with Hq = 70 km s~^ Mpc~^, O^ = 0.3, and Qa = 0.7. 



2. The SPIDER sample 



2.1. SDSS ugriz sample of ETGs 

The sample of galaxies is selected from the Sloan Digital Sky Survey (SDSS) DR6 in the 
redshift range 0.05 to 0.095 and with Mr < —20, where Mj. is the k-corrected S DSS Petrosian 



magn itude in the r-band. The k-correction is estimated using the software kcorrect ([Blanton et al 
20031 ). through a rest frame r-band filter blue-shifted by a factor (1 + zq), where zq = 0.0725, the 
median redshift of the ETG sample (see Paper I , Section 3.1). The lower redshift limit is chosen to 



minimize the aperture bias (IGomez et al.ll2003l). while t he upper redshift limit guarantees a high 
level of completeness (according to lSorrentino et al.ll2006l ) and allows us to define a volume-complete 
sample of bright early- type systems. ETGs foUow t wo different trends in the size- luminosity diagram 



(jCapaccioli et al. 



1992 



Graham Sz GuzmanI l2003l ) . The separation between these two families of 
bright and ordinary ellipticals occurs at an absolute -B-band magnitude of —19, corresponding to the 
magnitude limit of M^ ~ —20 adopted for this selection. At the upper redshift limit of z = 0.095, 
the magnitude cut of —20 also corresponds approximately to the magnitude limit where the SDSS 
spectroscopy is co mplete (i.e. P e trosia n magnitude of rrir ~ 17.8), making the sample volume 
limited. Following iBernardi et al.l ()2003l ). we define ETGs using the SDSS spectroscopic parameter 
eClass, which indicates the spectral type of a galaxy on the basis of a principal component analysis, 
and the SDSS photometric parameter fracDeVr, which measures the fraction of galaxy light that is 
better fitted by a de Vaucouleurs (rather than an exponential) law. In this contribution, ETGs are 
those systems with eClass < and fracDeVr > 0.8. The SDSS selection criteria and completeness 
of the ETG sample, part of the SPIDER project, are further detailed in Paper I. 



2.2. UKIDSS YJHK photometry 

The SPIDER sample consists of 39,993 ETGs, with available ugriz photometry and spec- 
troscopy. This SDSS sample is then matched to the UKIRT Infrared Deep Sky Survey (UKIDSS) 
Large Area Survey DR4. The UKIDSS-LAS DR4 provides NIR photometry in the YJHK-hands 
over ~1000 square degrees on the sky, with significant overlap with SDSS. The YHK-hand data 
have a pixel scale of 0.4" pixel"^, matching almost exactly the resolution of the SDSS frames 
(0.396" pixel"^). J-band observations are carried out with a resolution of 0.4" pixel"^, and then 
interleaved to a subpixel grid, resulting in stacked frames with a resolution of 0.2" pixel"^. The 
yji^ i^-stacked images (multiframes) have average depths of 20.2, 19.6, 18.8, and 18.2 magnitudes 
(in the Vega system), respectively. For each ETG in the SDSS sample, we searched for the nearest 
UKIDSS detection within a radius of 1", considering only UKIDSS frames with good quality flags 
(ppErrBits < 16). Of these galaxies, 5,080 objects also have available photometry in the YJHK 
wavebands from UKIDSS. Hereafter, we refer to these samples as the complete and optical+NIR 
samples of ETGs, respectively. 



As detailed in Pape r I, all grizYJHK frames have been homogeneously processed with 2DPH0T 



(JLa Barbera et al.ll2008l ). an automated software environment that performs several tasks, such as 
catalog extraction (using SExtractor), star/galaxy separation, and galaxy surface photometry. For 
each ETG, magnitudes are measured within the same aperture in all ugrizYJHK wavebands. 
Unless otherwise stated, in the present work we use the ugrizYJHK magnitudes measured within 
an adaptive circular aperture of 3 x r^^i, where r^^i is the Kron radius in the i-band. We use 
i-band Kron radii because of the larger S/N ratio of i relative to zYJHK frames (see Section [5]), 
and its lower sensitivity to young stellar populations in a galaxy relative to ugr. All magnitudes 
herein are in the AB system. For more than 95% of all ETGs in the SPIDER sample, the 3 x r^,? 
aperture is at least three times larger than the seeing FWHM of the grizYJHK frames, making 
the Kron magnitudes essentially independent of the seeing variation from g through K (see Paper 
I). We also note that the first four papers in this series do not use u-band photometry, as the S/N 
of SDSS ti-band frames is too low to measure reliable s urface photometry. For t he present study, 



we have processed the u-band images with SExtractor (JBertin &: ArnoutsI 119961 ). using the same 
setup as for the grizYJHK wavebands, now resulting in 5,068 objects in the optical+NIR sample. 
All magnitudes have been corrected for Galactic extinction, as detailed in Paper I. 



2.3. SDSS spectroscopy 

SDSS spectra in the ran ge 3800— 9200A are analyzed with the spectral fitting code STARLIGHT 



(jCid Fernandes et al.l |2005| ) to determine galaxy ages, metallicities, and interstellar extinctions, 
among other parameters. STARLIGHT finds the combination of SSP models that, normalized 
and broadened with a given sigma, best matches the observed spectrum. To this eff'ect, objects 
are de-redshifted and corrected for foreground extinction, following the recipes described in Paper 
I. For the galaxy spectra in the SPIDER sample, the median value of the resolution varies from 
~2.8A (FWHM) in the blue (4000A) up to ~3.7A (FWHM) in the red (8OOOA). The theoretical 
spectra used for comparison (and described below in Section [3]) have similar resolution to that of 
SDSS spectra across the whole wavelength range from '~3000A to ~9000A and are therefore not 
resampled for this study. 



3. Stellar population models 

We begin by constructing a library of model spectra (SEDs), generated through different stellar 
population synthesis techniques that encompass a variety of stellar evolutionary tracks. The models 
span a wide range in age, metallicity, extinction, and IMF. In the following, we describe the SEDs 
used to fit the broadband colors as well as the spectra of ETGs. 



3.1. SEDs for fitting broadband colors 

We use ste l lar population synth esis models to convert galaxy luminosity into stellar mass (e.g. 



Bell et al. 



2003 



Fontana et al.ll2004l ). The stellar mass is derived from the factor needed to rescale 
the spectrum from the best fit theoretical stellar population (normalized at 1 Mq) to the intrinsic 
(observed) luminosities. The models we use to fit galaxy colors are described below. 



3.1.1. Wruzual & Charloi hOOA ) 



Among the variety of stellar evolution libraries provided by this code (hereafter BC03), we use 
models from the Padova 1994 library, which provide a median spectral resolution R = 2000 over 
the spectral range from 3200 to 9500A (STELIB), and R = 300 outside this ra- n ge (B aSeL 3.1; 
see BC03 for details). Models of stars with T > 50,000 K are taken from iRauchI (|2003l ). and the 
spectral models of the the rmally-pulsating asymptotic giant branch (TP-AGB) phase are based on 
Vassiliadis k. Wood! (|l993l ) . We utilize r nodels with fou r diff erent metallicitie s: 0.2Zq , 0.4:Zq , Zq , and 2.5Zq . 
Within this code, we also consider both lSalpeteiJ (|l955l ) and lChabrieiJ (|2003l ) IMFs, with mass limits 
of O.IM0 < M < lOOM©. Unless otherwise noted, we assume a Chabrier IMF due to its theoreti- 
cal motivation. The star-formation histories include exponentially declining SFR ex e~*' "^ with the 
values of r given below, no gas recycling, and a i = 20 Gyr cutoff time for star-formation. This 
ensures that no models have their SF cut off before an e-folding time. Table [U provides a summary 
of the stellar population model parameters, where the last column provides stellar ages used in 
comparing LePhare and STARLIGHT results (discussed in Section [5]). 



Table 1: Stellar population model parameters (BC03/CB10). 
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The SEDs were generated for a grid of 64 ages in the range 0.8—14.2 Gyr. Since the mixing of 
dust and stars in (early-type) galaxies is far from well understood, for the purposes of the present 



work we describe d ust extinction with a simplified approach, where attenuation is appHed to the 
templates using the ICardelli et alj (|l989l ) law, with E{B — V) in the range 0—0.5. The reddening 
E{B — V) is limited to 0.5 magnitudes to avoid incorrect fitting of observationally red galaxies 



as hig hly reddened blue galaxies; such objects should be absent from our sample. Illbert et al, 
(|2010l ) show that the stellar mass measurements are sensitive to the extinction law, with absolute 
median differences of 0.14 and 0.27 dex for two spectroscop i c sarn ples of high redshift galaxies 
when comparing the lCalzetti et al.l (|20od ) and lcharlot fc Falj (J200ol . CF2000) extinction laws. We 
present similar results in Section 17. 7[ Unless otherwise noted, we use the Cardelli extinction law 
throughout this paper. We note that most studies using SED fitting to determine stellar masses 
have used the Calzetti extinction law. However, this law is theoretically motivated for strong 
starbursts. While this may be reasonable for field surveys, where the majority of galaxies are likely 
to be star-forming, it is inappropriate for our sample, where we expect ETGs to be dominated by 
older stellar populations. Therefore, we utilize the Cardelli relation as our baseline. 



3.1.2. Chariot & Bruzual (2010) 



An updated version of this code, which includes the new prescription of iMarigo &: Girardi 



( 2007 ) for the TP-AGB evolution of low and intermediate-mass stars, is described in 



Bruzual 



20071 ) and is used in previous papers in this series (Papers II, III, IV). We consider its preliminary 



version in the following comparisons, referring to it as the CBIO code (Chariot &: Bruzual, private 
communication) . These templates use the STELIB library over the entire wavelength range covered 
by the ugrizYJHK bands. We use the same stellar population parameters as BC03, described 
above and listed in Table [TJ 



3.1.3. PEG A SE. 2: iFioc & Bocca- VolmeranaA 1(1991 ) 



The code by lFioc fc Rocca-Volmerangd (119971 ) is based on the same tracks and the same stellar 
spectral library as BC03. The difference between the spectrophotometric models is given by different 
stellar spectra assumptions for the hottest stars (T > 50, 000 A') and a different prescription for 
the TP-AGB phase , both affecting the r e sults at young ages (< 2 Gyr). Indeed, spectra of hot 
stars are taken from IClegg fc MiddlemassI (119871) , while the models of the TP-AGB phase are based 
on the prescriptions of iGroenewegen &: de Jong] ()1993l ). For comparison with the results obtained 
with the BC03/CB10 codes, we have used models with e xponentially declining SFR with time 
scales r = 0.1,0.3, 1.0, 2.0,3.0, 5.0, 10.0, and 15.0 Gyr and a IScalol (Il986l l IMF. Solar metallicity is 
assumed for each template, based on the readily available models within LePhare. The models 
are set up with no infall (i.e. all the gas available to form stars is assumed to be in place at time 
t = 0) and no galactic wind, and the default value of 0.05 has been assumed for the parameter 
representing the fraction of close binary systems. Because this library is available in the public 
version of LePhare, we choose to compare its results here. 



The SEDs were generated for a grid of 39 ages (in t he ra nge 0.8—14.2 Gyr). Dust extinction 
was apphed to the templates, using the ICardehi et al.l (|l989l ) law (again with E{B — V) in the 
range 0—0.5). 



3.2. SEDs for fitting spectra 

For each galaxy in the optical+NIR sample, we fit BC03/CB10 templates in the wavelength 
range 3800— 9200A to its SPS S spectrum. Per o ther studies using this sample, we also run 
STARLIGHT u sing SSPs from dVazdekis et al.ll20ia . hereafter, M09) that are based on the MILES 



stellar library (jSanchez-Blazquez et al. 



200 



3), 



which has an almost complete coverage of stellar 
atmospheric parameters, containing spectra of stars in the solar neighborhood. These SEDs cover 
the spectral range 3800— 7500A with a spectral resolution of 2.3A. Hence, they are well suited to 
analyzing SDSS spectra, whose spectral resolution is ~2.36A (FWHM) in the wavelength interval 
4800— 5350A. By improving on the instrumental homogeneity among stellar spectra, which share 
the same wavelength scale and resolution, and the accuracy of relative flux calibrations as compared 
to previous stellar population model catalogs, we expect M09 to contribute only small systematic 
biases in the estimates of spectroscopic parameters. Since support for the M09 library is not pro- 
vided within LePhare by default, and given the minimal usable wavelength range, we have opted to 
use BC03 models in the photometric SED fitting to compute the theoretical magnitudes, provided 
reliable spectroscopic measurements are obtained from the M09 models. 

To run STARLIGHT, we select a basis of 76 MILES SSPs, with ages and metallicities listed 
in Table [H All models hav e \a/Fe\= (i . e. sol ar abundance ratio). Dust extinction was applied 
to the templates, using the ICardelli et al.l (jl989l ) law. 



4. SED fitting using a x^-minimization 

The theoretical SEDs are compared with either the observed galaxy colors or the observed 
spectra to determine the age, metallicity, and color excess E{B — V) that best fits the observations. 
In the case of galaxy colors, the stellar population models (see Section [3|) can be used to produce 
theoretical S EDs at d i fferen t redshifts, and hence the fitting provides an estimate of the photometric 
redshift (see iBenited (J200d ) for a review) . 



4.1. Fitting broadband pliotometry with LePhare 



We perform the SED fitting with the photometric redshift code LePhare (jArnouts et al.lll999l : 
Ilbert et al.ll2006l ). which uses a ^^-minimization technique. Other po pular photometric redshift 
codes have been tested on various figures of merit. lAbdalla et al.l (|2008l ) use ~13,000 spectroscopic 



redshifts from luminous red galaxies in SDSS DR6 to test six publicly available p hotometric redshift 
codes, including LePhare and the popular codes HyperZ (JBolzonella et al.ll2000l ) and BPZ (jBenitez 
2000), finding that LePhare performs best in the lower redshift intervals. This result and the 
existing collaboration with the code developers motivates this choice for our study. Each SED 
is redshifted up to z = 0.3 in steps of 6z = 0.005 and convol ved with the SDSS/UKIDSS filter 
transmission curves Q The opacity of the intergalactic medium ( Madaulll995l ) is taken into account 
within LePhare. The merit function x^ is defined as 



x\z,T,A) 



E 



?f 



FL-AxFL,{^,T) 



f 



a 



f 

obs 



(1) 



predC^^ z) is the fiux predicted for a template T at redshift z. F^^^ 
^obs i^ ^^^ associated error, converted from the AB system. The index / refers to the specific filter 
and Nf is the number of filters. The photometric redshift is estimated from the minimization of x^ 
varying the three free parameters z, T, and the normalization factor A. This normalization factor 
depends on the choice of waveband used for scaling and is calculated as 



A 



E 



obs pred 

i^bs)' 



^f 



/=1 



pred 



(2) 



where / refers to the waveband(s) used for scaling. Studies often use near-infrared bands (e.g. 
iiT-band) for scaling, since it is only we akly affected by dust extinction and is quite i nsensitive to 



the p resence of young, luminous stars (JLilly &: Longaiiill984l : iGlazebrook et al.lll995l : ICowie et al 



19961 ) . In Section 17. 4^ we present a comparison of galaxy stellar masses for several choices of band 
scaling. Here, we note only that it produces a negligible difference in the mass and, unless otherwise 
noted, we utilize all available wavebands to scale the SED. 



4-. 1.1. Improving data and model matching with photometric zero-point offsets 



Th e SDSS/UKipSS pho tometric zero-points are uncertain at the few percent level (see lFukugita et al 

(j 19961 ): iHewett et al.l (J2006l ). respectively). To insure further accuracy of data and model fiux 



matching, it is important to provide reliable uncertainties in the zero-point magnitudes. To this 
end, we utilize the complete and optical -l-NIR samples with spectroscopic redshifts, including ob- 
jects in the range 14.0 < i < 16.4 (brighter 25% of ugrizYJHK galaxies). Using a x^-minimization 
(Equation [T|) at fixed redshift, we determine for each galaxy the corresponding best-fitting stellar 
population template. LePhare notes in each case -FJ^^, the observed fiux in the filter f. A x FL.^^ 



■^The SDSS filter curves were obtained from'http : //www . sdss . orR/drS/i nstrumeiits/imaRer /f ilters/index . htrnlj 
which include instrument efficiency. UKIDSS filter curves are taken from iHewett et al.l ( 20061 ) . All transmission is 
considered at an airmass of 1.3. 



10 



is the predicted flux derived from the best fit template and rescaled using the normalization factor 
A of Equation [21 For each filter /, the sum 



N, 



gal 



r = yAU>< Hrea - fL + «0 /-oi (3) 



/-^ W pred "^ obs ^ "^ } I " obs 

is minimized, leaving s^ as a free parameter. For random, normally distributed uncertainties in the 
flux measurement, the average deviation s^ should be zero. Instead, we observe some systematic 
differences, which are listed for the optical+NIR sample at fixed redshift in Table [2] (median, 
converted to magnitudes. In our data, these differences never exceed 0.042 mag (F-band) and have 
an average amplitude of 0.019 mag, using BC03. We see that these differences depend very weakly 
on the magnitude cut, listed as 14.0 < mi < x, adopted to select the sample and are also almost 
independent from the set of templates used in the fitting. The sizes of these systematic differences 
are comparable to the level at which the SDSS/UKIDSS photometric systems differ from a true 
AB system. We then proceed to correct the predicted apparent magnitudes for these systematic 
differences, using the correction factor s^ . If we repeat a second time the procedure of template 
fitting after having adjusted the zero-points, the best fit templates may change. LePhare checks 
that the process is converging: after two iterations each estimated correction s-' varies less than 
0.0013 mag. Since the uncertainties in these zero-point corrections are not less than 0.01 mag, this 
error is added in quadrature to the apparent magnitude errors. 

Table 2: Systematic differences s^ between the observed and predicted fiuxes. 
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4.2. Spectral fitting 



For each spectrum in the optical-|-NIR sample, we run STARLIGHT (|Cid Fernandes et al 



20051 ) using a suite of models with stellar population parameters given in Table [H separately for 
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BC03/CB10 and M09 (see Section [3.2p . The code uses a x^-minimization, with a figure of merit 
given by 

x\x,Av,v,,a,) = Y,[{Ox-Tx)u:x]\ (4) 

A 

where 0\ is the observed spectrum at wavelength A, T\ is the model spectrum, and (jJ\ is a weighting 
factor. The model spectrum is a convolution of the intrinsic spectrum x with a Gaussian filter 
centered at v^ with dispersion a^. Since both BC03/CB10 and MOO have similar spectral resolutions 
to the SDSS spectra, the a^ provides a direct estimate of the galaxy velocity dispersion (see Paper 
I for details). The intrinsic spectrum aj is a linear combination of a set of stellar population 
models (the base) provided as input to the code. The coefficients of this linear combination are 
computed from STARLIGHT as part of the x^-n^inimization procedure. For our purposes, we 
do not use this feature of STARLIGHT, but rather use the reduced-x^ statistic provided by the 
code for each model in the base to compute the age, [Fe/H], and Ay of the best fit (i.e. lowest 
X^) stellar population model. For processing efficiency the code permits a maximum of 300 base 
models, so to allow a reasonably dense grid of ages/metallicities (columns 3,4 of Tabled]), we chose 
to use single {t = 0) burst models, which is justified by the mostly short-burst fits shown in Figure 
[U (see top panels). Here we plot the distribution of best fit BC03/CB10/PEGASE.2 models for 
the optical+NIR sample, where the abscissa shows increasing r (i.e. decay time scale for star- 
formation), with every fourth value being a constant metallicity - for reference, the first five peaks 
(BC03/CB10) represent solar metallicity models with r < 3 Gyr. A comparison of best fit ages, 
metallicities, and extinctions between LePhare and STARLIGHT is given in Section F5. 21 



5. SED fitting reliability 

In this section, we compare results obtained by fitting broadband colors and spectra of ETGs. 
We begin by comparing photometric and spectroscopic redshifts in Section [5. 11 while in Section [5^2 
we compare age, metallicity, and internal extinction estimates between LePhare and STARLIGHT. 



5.1. Accuracy of photometric redshifts 

We compare spectroscopic and photometric redshifts for the optical+NIR sample in Figure 
[2j The galaxies are binned to include 200 objects per 5zspec bin, and the error bars denote the 
dispersion defined by 1.48 x median(| Azj/(1 + Zspec))- This measurei nent of the scatte r corresponds 



to the rms for a Gaussian distribution and is unaffected by outliers (jllbert et al.ll2006l ). We recover 
~73% of the galaxies (BC03) with outlier rate, defined here as rj = \Az\ < 0.025(1 + Zgpec)- The 
accuracy of t hese photoinetric r edshifts improved with addition of the u-band data by ~7% in the 



outlier rate. Illbert et al.l (|2006l ) use their sample of '--'3,000 galaxies to confirm the importance of 
the u-band for low redshift galaxies, recovering ~80% of the photometric redshifts at z < 0.4 and 
~95% using the u-hand. Figure [2] shows the binned redshifts, using the Chabrier (Scalo) IMF for 
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BC03/CB10 (PEGASE.2) stellar population models, and the Cardelli extinction law. Photometric 
redshifts are well recovered, given the small redshift range, with outlier rates of ~73%, ~75%, ~76% 
for BC03/CB10/PEGASE.2, respectively. However, at higher redshifts, CBIO models produce 
systema t ically lower photometric redshifts. This is surprising given the improvements described in 



Bruzuall (|2007l ). although it does not change the results that follow, which use the spectroscopic 
redshift. We note that the extremely small redshift range of our sample makes these photo-z tests 
unsuitable for drawing conclusions for typical imaging surveys. 



5.2. Comparing photometric & spectroscopic measurements 

To compare photometrically and spectroscopically determined results, we use the same model 
parameters between LePhare and STARLIGHT as much as possible. All theoretical parameters are 
consistent between the two packages, except for the interstellar extinction, which is calculated on- 
the-fly within STARLIGHT. The comparison assumes that there is essentially no stellar population 
gradient between the spectroscopic (fiber) and photometric (Kron) apertures. This assumption is 
motivated by the lack of color gradients between fiber and Kron apertures, as shown in Appendix lAl 
As shown in Figure El there is good agreement between the distributions of photometric and spec- 
troscopic ages using CBIO (see bottom-left panel), with a median difference of 0.75 Gyr, whereas the 
sample median age (measured from AGE_MED in LePhare) |j estimated from the BC03 templates 
is 1 Gyr too low, slightly improved by addition of the w-band data. However, the new treatmen t 



of the TP-AGB phase yields isochrones up to 1 magnitude brighter in the iiT-band ([Bruzuall 120071 ) . 
where the dominant flux of these galaxies is expected - this would produce better fits between the 
photometric SED and older stellar populations. Despite this discrepancy, ugrizYJHK photometry 
reproduces the age-metallicity distribution much more accurately than the ugriz-hands alone, as 
expected, and so is used in the tests that follow. Both STARLIGHT and LePhare produce similar 
age/metallicity/extinction distributions for our optical+NIR sample, despite large differences for 
individual galaxies and strong systematic trends as a function of STARLIGHT output parameters 
(see Figured]). 



6. Method to estimate galaxy stellar masses 

In most cases, all 9 bands were used to compute the stellar mass. We provide the rescaled 
template stellar mass (measured from MASS_MED in LePhare) derived from 

logX = log{M,/Cx) + OAkcorx + 2logdpc - 2.0 + 0.4Mr" - 0.4mA', (5) 



^Stellar ages and masses are measured from the median of the hkeUhood function exp(— x^/2) rather than the 
best x^- This reduces stochastic mass errors due to individual models which happen to have very low x^ values, since 
there remains some degeneracy in the SED fits even with 9 filters. We also find that this method produces masses 
that are more consistent with other studies. 
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where J^^j L is the stellar mass-to-light ratio in the chosen filter centered at A, kcor\ is the k- 
correction, dpc is the cosmology dependent luminos ity distance in pc, and mv is the apparent 



magnitude in the observed filter centered at A' (see iLonghetti &: Saraccd (|2009l ) for a derivation 
of Equation [5]) . As with stellar age, this measurement is taken from the median of the likelihood 
function exp(— x^/2). Note that stellar masses computed in this manner are only partial due to 
the adaptive aperture of 3 x rK,i\ total stellar masses would require an aperture correction, which 
is irrelevant for this type of study (however, see Section [8]). We then compute the median and its 
asymmetric la uncertainties for histograms of A (log }A) for differences resulting from the choices of 
stellar population model, initial mass function, interstellar extinction law, and photometric bands. 

We also determine stellar mass, using the optical+NIR sample, for fixed age, metallicity, 
and extinction, derived from STARLIGHT, in Sections 17. 7117. 81 We define the difference A as 
log A4 fixed — ^og Ai free for each parameter. Operationally, fixed refers to constraining the specified 
parameter and leaving all other parameters (except for redshift) free in the x^-fitting. This allows 
us to quantify the accuracy of a stellar mass estimate when these parameters are unknown, and 
then derive corrections for any offset. Therefore, we define a mass correction term, e, which can 
be applied to the stellar mass measurements over a range of photometrically determined ages and 
metallicities, for a given best fit template (i.e. through a linear regression line). This term is 
computed as the median of each A(logA^) histogram described above for fixed-parameter minus 
free-parameter. In Section [?.8. 21 we discuss the use of this parameter in correcting the stellar mass 
distribution of our samples. 

Lastly, we present two comparisons of our measured stellar masses to similar estimates from 
a group at the Max Planck Institute and to those presented in Paper IV ~ this is discussed in 
Section [HI We note here that the overall agreement with separate photometrically estimated stellar 
masses (MPI) is excellent, showing a negligible offset with 2a scatter of about 0.1 dex. We find 
that inclusion of short-burst model SEDs when fit with a young (<3 Gyr) stellar population yield 
stellar masses that are in most cases ~0.2 dex too low. 



7. Biases and uncertainties in stellar masses 

Unless otherwise noted, the results presented herein use the optical-|-NIR sample of 5,068 
ETGs and their spectroscopic redshifts. In Sections 17.11 and 17.21 we test the impact of sample 
photometric errors and redshifts, respectively, on the resulting stellar mass estimates. The tests 
presented in Sections 17.2117. 71 were conducted separately for both BC03 and CBIO stellar population 
models. In Sections I7.7H7.81 we apply results from the STARLIGHT spectral fits to constrain the 
X^-minimization, using the BC03, CBIO, and M09 libraries, separately for each test. When spec- 
troscopically constraining t, Z, or Ay with the M09 library, we use BC03 models within LePhare, 
as described in Section 13.21 In Table [H we also summarize the results of the following tests for 
BC03/CB10 stellar population models. Similar results using the complete sample of ETGs are 
provided in Table [H The reader should take results from this sample with caution, considering the 
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caveats introduced by use of only five optical bandpasses and the discussions in Section 15.21 and 
below. Further, we do not perform the SED fitting for fixed age/metallicity/extinction on the com- 
plete (ugriz-only) sample due to resulting catastrophic errors in the constrained x^-minimization. 



7.1. Dependence on sample magnitude errors 

The photometric measurements of the complete and optical+NIR samples in this study have 
high S/N ratios (see Figure ED . To investigate the dependence of galaxy stellar mass estimates 
on photometric errors, we used a Monte-Carlo approach to varying the ugrizYJHK magnitudes. 
Each magnitude was varied according to a normal distribution, centered on the observed apparent 
magnitude, using the la error limits provided by SExtractor. For fixed redshift, the resulting 
differences in stellar masses from the original sample (optical-1-NIR) are shown in Figure [6j The 
masses show a dispersion of ±0.12 dex at the 95% CL, which lie within uncertainties induced by 
the photometric redshift, but are surprisingly large given the S/N of this sample, with a negligible 
systematic bias introduced in the mass. 



7.2. Dependence on photometric redshift 

We quantified how the stellar mass accuracy is affected by the use of a photometric redshift, 
which is the typical case for a large imaging survey. Figure \7\ shows the difference between the 
stellar masses computed with the spectroscopic vs. photometric redshifts. The sample uses the 
entire redshift range 0.05 < z < 0.095, and all bands are used to scale the mass. We find a median 
difference of ~0.01 dex, with errors of +0.804/— 0.500 dex at the 95% CL, using BC03 models. 
This scatter is slightly larger than the systematic uncert ainties expected (0.2 dex at the 68% CL) 



in th e stellar mass estimate due to photometric redshifts (jPozzetti et al.ll2007l : iLonghetti fc Saracco 
20091 ) at low redshift Q For the optical+NIR sample, the average positive 2a limit on the photo-z 
is roughly 40% - at the largest redshift {z = 0.095, i.e. di = 435 Mpc), this yields a luminosity 
distance error of nearly 200 Mpc, or a stellar mass error (from Equation [5]) of ~0.4 dex. The 
photo-z also impacts the template best fit ages differently, depending on stellar population model. 
We concluded that the unusually low redshift of our sample results in abnormally large mass errors 
(compared to those reported in the literature) due to the use of photo-zs, since the fractional errors 
in the luminosity distances are large. For all remaining tests, we use the spectroscopic redshifts. 



Pozzetti et al.l (|2007f ) find uncertainties of ^-^0.2 dex at 2: < 0.4, with increasing errors at lower redshifts. Our 



redshift range is much lower, so these errors may be expected. 
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7.3. Dependence on stellar population model 

The comparison of stellar masses between different models is shown in Figures [8] and [9] (see top- 
left panels). To allow a meaningful comparison with BC03/CB10 and PEGASE.2 models, here we 
include only objects in the optical+NIR sample that are fit with a solar metallicity {Nbcos = 3, 371, 
NcBio = 4, 240) 1j. The stellar mass comparisons between models yield small offsets. CBIO models 
produce slightly higher masses by about 0.005 dex, with a scatter of ~0.14 dex with respect to 
BC03 masses using the spectroscopic redshift (see Table H]). This result seems to be in disagreement 
with, instead, BC03 models having systematically higher NIR A4^:/C than CBIO models, at given 
age and metallicity (see e.g. lEminian et al.ll200 8i). Indeed, LePhare gives lower photometric ages 



when using BC03 (see Figure [3]), resulting in lower M.^1 L ratios when compared to CBIO results 
This might compensate for any variation of A^*/£ due to the different treatment of the TP-AGB 



phase between BC03 and CBIO. This seems to hold only for ETGs with old ages, as lMuzzin et al 



(J2009l ) actually found CBIO (as of 2008) based stellar masses to be sy stematically smaller than 
those obtained from BC03 for high-redshift (z ~ 2.3) galaxies. Moreover, Ivan der Wei et al.l (J2006l ) 
reported that, for z ~ galaxies, including the TP-AGB phase, as done in iMarastonI (J2005l ) stellar 
population models, leads to lower stellar masses than those obtained with BC03, when fitting both 
optical and NIR broad-band colors, more consistent with what is reported here. On the other 
hand, PEGASE.2 models produce systematically lower masses than BC03/CB10. The offsets in 
these comparisons are likely due to the different prescriptions for th e luminous stars dom i nating 
the iC-band flux, and are comparable to other factors discussed below. iLonghetti &: Saraccd (|2009l ) 
measure the effect of different TP-AGB prescriptions (assuming these stars dominate the -fC-band 
flux) between stellar population models on the resulting stellar mass when the age is known, finding 
agreement within at least 20%, between BC03, CBIO, and PEGASE.2 models. 



Earlier studies (Pozzetti et al. 



2007 



Treu et al.ll2010l ) find that a Chabrier IMF underestimates 



galaxy stellar masses b y nearly a factor of 2, compared to those derived with the Salpeter IMF. 
Cappellari et al.l (|2006l ) used different models than these authors and find that galaxy stellar masses 
based on a Salpeter IMF were in some cases too high compared to those determined with stellar 
kinematics, reaching the same conclusion. We use the stellar population models of BC03/CB10 to 
further compare model dependent stellar masses (with star- formation histories, metallicities, and 
color excesses given in Tabled]), using different IMFs. Figures [8] and [9] show the dispersion in 
A(log A^) between the Chabrier and Salpeter IMFs (see top-right panels). The zero-point offset is 
—0.227 (—0.225) dex using BC03 (CBIO) models, or a factor of ^2^ which differs from the analytical 
estimate (see Appendix [B]) by ~0.1 dex. 



^These results assume a lScalj (llOSd ) IMF for PEGASE.2 models. To compare with BC03/CB10 results, which 
use a Chabrier IMF, we remove the bias in stellar mass introduced by the IMF, as described in the Appendix iBl 
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7.4. Dependence on photometric bandpasses fitted 

Here we highlight the systematic effects of band scahng in the x^-niinimization as well as inclu- 
sion of ultraviolet and NIR continua, separately, in the SED fitting. Again, we remind the reader 
that these results apply specifically to a low redshift population of ETGs, where the continuum 
levels behave differently with redshift for ETGs and especially for other morphological types. 

In Figure [TUl we plot the differences in galaxy stellar mass for several choices of band scaling 
- ugrizYJHK, rK, and HK. For our optical+NIR sample of low redshift ETGs, this allows for 
a measurement of the sensitivity of the overall luminosity scaling of the SED to magnitude errors, 
specifically targeting the K-h&nd. The median, dispersion, kurtosis, and skewness for these stellar 
mass differences are given in Table [3l suggesting that there is a systematic bias in the stellar mass 
of 0.005 dex when using all bands as opposed to using only the H and K bands. The differences 
are sharply peaked around zero, as indicated by the high excess kurtosis. 

Table 3: Characteristics of A(log Al) for ugrizYJHK — xK (with x = r, H) band scaling. 



Median32CT Kurtosis Skewness 



ugrizYJHK -rK -8.000E-3[J:K 9.09 -0.734 

ugrizYJHK -HK -4.500E-3|];K 5.61 1.00 



These offsets are small compared to errors introduced by other factors explored later in this 
paper. The excess of positive residuals for ugrizYJHK — HK scaling is possibly due to a known 
poor modeling effect in the near-infrared described in Sections and 17.31 which predicts if -band 
magnitudes that are too faint. Therefore, we use all bands to measure the stellar mass, so that the 
scale factor is less sensitive to errors in any one given band. In the minimization, the best x^ ^^ 
each redshift step is saved to build the function F{z) = ex p[— y^^y, (z) /2] . This function is used to 



refine the photo-z solution with a parabolic interpolation ()Bevingtonl Il969l ) . 



7.5. Dependence on inclusion of rest-frame UV photometry 

We use the optical-|-NIR sample to compare stellar masses computed with and without u-band 
data, for BC03/CB10 models. We find scatters of ~0.14 and ~0.12 dex, respectively, in A(log A^), 
as shown in Figures El and El with negligible systematic bias in the mass (see green histograms in 



bottom-right panels). IConrov Sz Gunnl (J2010l ) describe the infiuence of the 4000A break (quantified 
as D„4000) on ugr colors for red sequence galaxies and find that BC03 models, among others, 
predict D„4000 strengths that are generally too large. It is possible that this effect, which alone 
would yield older stellar ages in the x^-fitting, competes with the underpredicted NIR luminosities 
to produce the scatter observed in Figure 13 (for u—no-u). This uncertainty would be difficult 
to quantify, as these models suffer in reproducing the age-metallicity plane found from spectral 
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fits. However, with the modifications in CBIO models affecting spectra redward of ~6000A, the 
u—no-u stellar mass differences in Figure [9] have a lower systematic bias and dispersion than the 
BC03 models, showing that CBIO treats optical+NIR spectral regions more consistently with the 
observations. Hence, including u-band data simply provides more robust stellar mass estimates. 



7.6. Dependence on rest-frame NIR photometry 

We also test the absence of near-infrared data on the stellar mass estimate (see bottom-right 
panels of Figures [S] and E|) • The A(log A^) offset in ugrizY JHK — ugriz stellar masses is approx- 
imately -0.057 ± 0.10 (-0.042 ± 0.12) dex for BC03 (CBIO) models, with n^riz-only photometry 
producing higher stellar masses in both cases. iPozzetti et al.l ()2007l ) find a similar result - that 
is, higher masses by ~0.1 dex using optical-only photometry - for their i^-selected spectroscopic 
sample of high redshift galaxies. Given the changes in TP-AGB evolution between BC03/CB10 
models, these relatively smaller dispersions suggest that the 4000 A break is a more sensitive con- 
straint in the SED fits than is the continuum level of the NIR data for this sample of low redshift 
ETGs. 



7.7. Dependence on interstellar extinction 



We compare stellar masses derived using the ICardelli et al.l (119891 ) and ICalzetti et al.l (120001 ) 
extinction laws with E{B — V) values listed in Table [TJ These empirical laws are line-of-sight 
dependent with average Ry equal to 3.1 and 4.05 ± 0.80, respectively, in popular photometric 
redshift codes (e.g. BPZ, Hyper Z, LePhare). Using BC03 models, this mass difference yields errors 
of +0.07/— 0.14 dex for the optical+NIR sample with negligible systematic bias in the mass. These 
extinction laws produce excellent mass agreement, however, using CBIO models, which might be 
related to the fact that CBIO models provide SEDs that are more consisten t with the obser v ation s 
in both the optical and NIR. We also measure stellar masses using the ICharlot &: Falll (120001 ) 
extinction law, which is provided with the GALAXEV code and requires as inputs the total 
effective F-band optical depth Ty and the fraction /i of it contributed by the interstellar medium. 
The default values of Ty = 1.0 and // = 0.3 are assumed for this test, where it should be noted that 
these values best describe a population of star-forming galaxies. Our optical+NIR sample yields 
a A(logA^) zero-point offset of —0.10 dex for BC03 models, with CBIO models producing slightly 
better agreement. Note that this difference is comparable to the offsets found bv lllbert et al.l (J2010l ) 
for the z ~ 1 zCOSMOS sample (A(log Al)(Caizetti-CF2000) is —0.14 dex), where they also show 
that this systematic offset is larger for massive galaxies with a high SFR. 



Finally, we utilize spectroscopic measurements of Ay (from STARLIGHT) and the lCardelli et al 



(jl989l ) extinction law to constrain the interstellar extinction for each galaxy. In LePhare, we fix 
the color excess E{B — V) based on the relationship E{B — V) = Ay/Ry. We find that for our 
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optical+NIR sample, the uncertainties in the stellar mass estimate due to unknown extinction 
are comparable (~0.2 dex for BC03) to the effect of using different extinction laws. However, 
M09/CB10 models yield much larger uncertainties of ~0.4 dex. We find a negligible systematic 
trend in mass difference as a function of spectroscopically determined extinction for all models. 
Figure [11] shows however that when fixing the extinction to the spectroscopic value measured from 
M09/CB10 models, stellar masses are systematically lower than the free parameter case by ~0.15 
dex for all Ay values. A linear regression fit to the trend of A(log7W) offset with amount of 
extinction (in magnitudes) produces 

(0.132 ±0.832)Ay + (-0.156 ±0.263) for M09, 

A(log7W) = I (-6.50E-3 ± 0A25)Av + (8.54E-3 ± 0.140) for BC03, (6) 

(0.0471 ± 1.30)^y + (-0.132 ± 0.330) for CBIO 

where the la uncertainties on the regression coefficients are provided above. Large uncertainties in 
A(log A^) affect the reliability of this fit, but it is provided for reference. 

Fitting broadband photometry seems to produce lower interstellar extinction than spectral 
fitting for our sample of ETGs, as shown in the right panels in Figure [H This would lower the 
0.4 X my term in Equation [5] and produce a higher stellar mass estimate. We note here that we 
do not provide a mass correction analysis as a function of photometrically determined extinction, 
due to the shallow trend in A(log A^) with Ay- Furthermore, the Ay (or equivalently E{B — V)) 
parameter space cannot be specified directly within the STARLIGHT x^-fitting. 



7.8. Dependence on stellar age and metallicity 

This section spotlights both the effect of different best fit ages and metallicities on the resulting 
stellar masses, and suggests how blue continuum fiuxes and mass-to- light ratios could explain the 
large stellar mass differences due to photometrically determined ages. To this end, we have used 
optical spectra to constrain independently age and metallicity in the models. 



7.8.1. Dependence on stellar metallicity 

The Padova 1994 stellar evolutionary tracks used by BC03/CB10 encompass six metallicities 
in the range Z/Zq = 0.005 — 2.5, with the four most metal-rich models chosen for this study. The 
M09 model fits use Z/Zq = 0.2, 0.4, 1.0, 1.6 metallicities from the Padova 2000 library. We separate 
galaxies on the basis of the metallicity ranges given in Tabled] and run LePhare with the metallicity 
constrained to the closest value. The result is a systematic offset for unknown metallicity consistent 
with zero for most objects fit with M09 models. However, metallicity is more poorly constrained 
with BC03/CB10 models, yielding a larger systematic offset in the stellar mass, with uncertainties 
as high as ~0.2 dex at low metallicity (see Table Hj). The linear regression fit to our optical+NIR 
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sample yields the following functional form: 
A(log7W) = 



(-0.138 ± 0.248)Z + (0.0758 ± 0.136) for M09, 

(-0.162 ±0.147)Z + (0.0241 ±0.0769) for BC03, (7) 

(-0.320 ± 0.206)Z + (0.0296 ± 0.0362) for CBIO 



where the la micertainties on the regression coefficients are provided above. Large uncertainties 
in A(log A^) for the non-solar metallicity models affect the reliability of this fit, but it is provided 
for reference. 

We can see how the distribution of stellar masses changes for our sample by applying an offset 
to each mass as a function of its photometrically determined metallicity, which we call the mass 
correction, e. For galaxies with a given photometrically-determined metallicity, we add an offset 
calculated from the linear fit in Figure [12] to the stellar mass measured with only the spectroscopic 
redshift constrained. Figure [T3l shows the distribution of stellar masses for our optical-!- NIR sample, 
corrected for these offsets. This technique shows that the median stellar mass changes by ~0.1 dex 
when using the M09/BC03 models. The distribution remains more consistent in CBIO models, 
which may be expected given the slightly better agreement between spectroscopic and photometric 
metallicities using these templates (see CBIO/Metallicity plot in Figured]). 



7.8.2. Dependence on stellar age 

As described in Section [7^ stellar masses of low redshift ETGs have uncertainties of ~0.14 
dex at the 95% CL due to the stellar population models alone. Furthermore, the photometrically 
determined stellar age is model-dependent. The relationship between dominant stellar population 
age and stellar mass is not easily quantifiable, but errors in the age are propagated to errors in the 
mass estimate. To this end, we have used stellar age measurements from STARLIGHT spectral fits 
to fix the age of the template used in the x^-minimization. 

Figure [TT] shows an estimate of the systematic bias in stellar mass for unconstrained age plotted 
as a function of spectroscopically determined age. We note that when constraining the stellar age in 
LePhare to its spectroscopically determined value, only 1,314 (1,810) galaxies, using BC03 (CBIO) 
models, have non-zero exp(— x^/2) values. There is a clear trend of increasing mass difference 
with age. Furthermore, using BC03 models in the SED fitting (provided spectroscopic constraints 
separately from M09/BC03) yields systematically lower stellar masses at all ages. This trend is 
also evident from CBIO models, though stellar masses for unconstrained age are systematically 
higher below about 4 Gyr. The linear regression fit to the optical-!- NIR sample yields the following 
functional form: 



A(log7W) 



(0.0359 ± 0.0126)t + (0.0175 ± 0.0963) for M09, 

(0.0395 ± 0.0110)t + (-0.0148 ± 0.0837) for BC03, (8) 

(0.0467 ± 0.0237)t + (-0.144 ± 0.133) for CBIO 
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where the la uncertainties on the regression coefficients are provided above. Notice that there 
is a significant positive trend in stehar mass underestimation with increasing age that is weU 
approximated by the hnear model. 

This gross underestimate in mass for old stellar populations might be due to the fact that 
models underestimate the ages - possibly because of age-metallicity-extinction degeneracy or a 
failure of the models to reproduce the observed spectra - and then the SED fitting code increases 
the luminosity scale factor accordingly. To illustrate this, we select a single galaxy with spectro- 
scopically and photometrically determined ages of 12.5 Gyr and 2.2 Gyr, respectively, and compare 
the best fit redshifted, reddened photometric SEDs for free and fixed age (see Figure [T^. where the 
fixed age stellar mass is higher by -^0.42 dex. The galaxy is best fit by a Z/Zq = 0.4 metallicity, 
and differs by 0.1 in E(B — V), using BC03 models. Following Equation [5l the median difference 
in 0.4 X mx' in the g-hand is 0.06 mag (with higher flux in the 2.2 Gyr SED) and indistinguish- 
able from our broadband photometry alone. Older and younger populations have very different 
stellar mass-to-light ratios, where for instance the difference in log(A4*/£B) for the two best fit 
ages here is ~0.67 dex (with lower i?-band mass-to- light ratio in the 2.2 Gyr SED). As expected, 
the 4000A break is the most significantly differing feature in the broadband colors of these SEDs, 
with A\log{A4^/ Cb)] — A(0.4 x nig) accounting for most of the measured difference in stellar mass 
for this particular object, yet it cannot be photometrically distinguished despite the high S/N of 
our objects at this redshift |j. Furthermore, including the i^-band photometry seems to bias most 
SEDs towards young ages, assuming our ellipticals are indeed dominated by older stars. So, we can 
attribute larger differences in the stellar mass due to unknown age to the blue A1*/£ for different 
stellar population ages. It is also true that differing spectroscopic and photometric metallicities 
will produce different stellar masses, but given the trends provided in Section [7.8.11 it is likely that 
large differences in age dominate the effect of age-metallicity degeneracy on the galaxy stellar mass. 

We can see how the distribution of stellar masses changes for our sample by applying an offset 
to each mass as a function of its photometric age, similar to the procedure described at the end 
of Section 17.8.11 Figure [13] shows the distribution of stellar masses for our optical -l-NIR sample, 
corrected for these offsets. This technique shows that the median stellar mass changes by as much 
as ~0.1 dex when using the updated CBIO models. 



8. External comparison of stellar masses 

As a comparison to related work, we plot stellar masses obtained from our SED fits against 
those derived from a group at the Max Planck Institute (see Figure [T5]) . using the complete sample 



®We compare fluxes and mass-to-flght ratios in g and B-bands, respectively, for tliis object because M,/£ is 
provided for B and V^-bands from the BC03 code. Sloan g and Johnson-Cousins B fllters span roughly the same 
wavelength range and are just redward of the 4000A break, which due to the large difference in best fit stellar age, 
most constrains the continuum level of the model spectrum. 
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and BC03 stellar population models |3. We convert to total stellar mass by dividing each set of 
measured stellar masses by an aperture correction, corresponding to the fraction of z-band light 
contained within the z-band Kron aperture {LePhare) and the SDSS fiber (MPT). This procedure 
assumes that the stellar mass in a galaxy is distributed in the same way as the z-band light. We 
notice two trends in this comparison, with ~75% of the galaxies lying above an arbitrary line drawn 
0.1 dex below the 1:1 line. The subset of galaxies that lie above this line are in excellent agreement 
with those from MPI, with a negligible offset and 2a uncertainty of ~0.1 dex, likely due in part 
to the difference between how Kron and SDSS FiberMags are measured (e.g. correcting the latter 
for emission lines) as well as potential differences in model parameters. The distributions of best 
fit age, star-formation decay time scale, metallicity, and extinction are plotted in Figure [16] for 
galaxies lying above {top panels) and below (bottom panels) the dotted line drawn in Figure [T5l We 
see that with metallicity and extinction being similar for the two samples, the contribution from a 
high t/r is a redder galaxy - that is, such a galaxy has enough e-folding times to allow the younger 
stars to fade. This produces a degeneracy in the shape of the SEDs, as compared to a similar t/r 
for an older galaxy, resulting in a lower mass for the younger SED. We have directly verified this 
degeneracy for several objects with similar t/r, metallicity, extinction, and redshift, but different 
ages. 

We also compared the stellar masses that we present here with those presented in Paper IV, 
finding trends similar to those described above. Specifically, we find a negligible offset (0.003 dex) 
with a modest 2cr scatter of ~0.1 dex for galaxies (65%) lying above a line, drawn similarly to 
the dotted line in Figure [151 The stellar population model differences between the stellar masses 
presented in this paper ~ i.e. for BC03 models, with a Chabrier IMF, Cardelli extinction law, and 
spectroscopic redshift - and those presented in Paper IV is the permitted model decay time scales 
for star-formation (r), where in this paper we also include r = 0.1, 0.3, and 30.0 Gyr models, and 
use of the Calzetti law in Paper IV (which induces negligible offset in the stellar mass). Another 
difference is that in this paper we use all available bands to scale the SED, so that the scale factor 
is less sensitive to errors in any one given band. In Paper IV, we used only i^-band to scale the 
SED, due to its supposed independence from the effects of dust extinction and young, luminous 
stars (see Section HT]) . 



9. Summary & conclusions 

We measured systematic errors in galaxy stellar masses due to different ingredients in widely 
used models from stellar population synthesis techniques. Using a large sample of high S/N, low 
redshift ETGs, we want to accomplish two primary goals 



^However, detailed parameters, such as IMF, star formation history, extinction law, etc. are unspecified on their 
website (see |http : //www . mpa-garching . mpg . de/SDSS/DR7/mass_comp . html | , with differences likely contributing to 
the observed scatter. 
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- Determine at a high-level of significance, the statistical biases and uncertainties inherent in 
a stellar mass estimate; 

- Provide robust corrections to a sample of photometrically determined galaxy stellar masses, 
when spectroscopy is unavailable. 

We first presented a comparison of BC03, CBIO, and PEGASE.2 stellar mass estimates, finding 
that masses have dispersions of ~0.10 — 0.15 dex at the 95% CL, with this uncertainty arguably 
being due to the different treatment of the TP-AGB stellar evolutionary phase. In general, our tests 
find no statistically significant systematic mass bias at the 2a level due to variations of the stellar 
population models, extinction laws, and photometric bandpasses. This includes varying scaling 
bandpasses in the fitting technique as well as the magnitudes used in generating the photometric 
SED. However, a systematic offset of about —0.23 dex (inconsistent with zero at the 2(T level) is 
observed betwe en stellar masses p r edict ed from Chabrier and Salpeter IMFs. The newer TP-AGB 



calculations bv lMarigo &: Girardil (J2007I ) affecting i^T-band magnitudes also improve the quality of 
It-band fits, where the dispersion in galaxy stellar mass measured with and without n-band data 
is improved by nearly 20% from BC03 to CBIO models. It is interesting to note that stellar mass 
estimates are more consistent with no near-infrared photometry than they are without u-band data, 
at least using BC03 models with this sample of ETGs. 

The observed systematic stellar mass biases and dispersions for our optical+NIR and complete 
sample of ETGs are summarized in Table [H As expected for this redshift range, using th e photo-z 



produces the largest systematic uncertainties in the stellar m a ss est imate. We find that a I Chabrier 



(|2003l ) IMF produces lower stellar masses than a ISalpeteiJ (|l955l ) IMF by about 0.227 (0.225) 
dex for BC03 (CBIO) models, in agreement with other studies, but lower by about 0.1 dex than 
the difference predicted by direct integration of the formulae over the mass range 0.1—100 M.q 
(Appendix [Bj) . For the optical+NIR sample, this difference falls within the 3cr confidence limits 
of the observational estimate (i.e. 0.23 +0.15/-0.10 for BC03), while for the complete sample, it 
is within the 4cr confidence limits for BC03 and highly inconsistent with zero (> Act deviation) for 
CBIO. 

An important part of this work is the investigation of systematic effects on the photometrically 
determined stellar mass when spectroscopic constraints on age, metallicity, and extinction are un- 
available, and to provide statistically robust stellar mass corrections as a function of these photomet- 
ric parameters for low redshift early-type galaxies. We used 5,068 objects with ugrizYJHK pho- 
tometry and spectroscopy to achieve this goal. Offsets from extinction constraints are the smallest 
of these systematic effects; this implies that the SED fitting is less sensitive to the coarse extinction 
grid employed by stellar population models when estimating the stellar mass. We find a negligible 
trend with Ay in stellar mass difference for the optical+NIR sample. However, there is about a 0.15 
dex bias towards higher masses for all extinctions when Ay is unknown, using M09/CB10 models. 
The dispersion in these mass differences is about ±0.5 dex (consistent with zero), slightly lower for 
BC03 models, yet larger than uncertainties introduced by unknown age/metallicity, mainly because 
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of the smaller number of galaxies in a given extinction bin. Fitting broadband photometry seems to 
produce lower interstellar extinction than spectral fitting for our sample of ETGs, as shown in the 
right panels in Figure [H This would lower the 0.4 x m^' term in Equation [5] and produce a higher 
stellar mass estimate. We comment that for general surv eys, it might be more appropriate to fit 



an SED with one extinction law (e.g. ICalzetti et al.1120001 ). and then refit with a more appropriate 



law once an approximate age is established. 

The positive trend in A(logA4) with increasing metallicity is simply a result of the inverse 
trend shown in the middle panels of Figure [H That is, spectral fits on average produce lower 
metallicities than broadband photometric SED fitting for these objects; and for any given age, a 
lower metallicity SSP has a lower mass. This is merely a reflection o f the failure of such a coarse 
metallicity grid to constrain the model in the SED fitting, as expected. iDressler fc ShectmanI (119871 ) 



show that the size of the 4000A break in spheroidal galaxies is quite insensitive to metallicity, so 
even provided a wider range of stellar population metallicities, our result should remain unchanged 
in this regard for (old) ETGs. We chose a larger grid of stellar ages for M09/BC03/CB10 models, 
ranging from 0.5 — 12.5 Gyr to constrain age in the x^-niini™ization. We measure remarkably 
lower stellar masses when age is unknown for all young/ old populations (except < 4 Gyr for CBIO 
models), with uncertainties of around 0.2 — 0.3 dex. Indeed, for our low redshift ETGs with older 
stellar populations, we would be underestimating the stellar mass content by a factor of 2 (for 
galaxies with t ~ 10 Gyr from BC03/CB10 models). We can argue that the increasing trend with 
stellar age in Figure [11] is directly due to the overall underestimate in photometrically determined 
age for all models (see left panels of FigureH]), whereas we argue in SectionE821 a difference of ~0.7 
dex in the best fit age can yield a stellar mass difference of ~0.4 dex, explained by varying M.^1 L 
ratios between different stellar age SEDs. Notice that unknown ages, resulting from a mismatch 
between the true star formation history of a galaxy and the simplified form used in SED fits (e.g. 
an exponentially declining SFR), can also affect significantly t he st ellar mass deter mination of 



high-redshift galaxies, as discussed, e.g., bv lPozzetti et al.l (J2007l ) and iLee et al.l (J2010l ) 



In conclusion, we find uncertainties in the galaxy stellar mass due to different stellar pop- 
ulation models, IMFs, and extinction laws that are much less than uncertainties introduced by 
using the photo-z over this small redshift range. More notable are offsets measured between 
Chabrier/Salpeter IMFs, Cardelli/CF2000 extinction laws, and changes in TP-AGB evolution- 
ary prescriptions - between CB10/PEGASE.2 models. The discrepancies in ages, metallicities, and 
extinctions between spectroscopic and photometric measurement techniques propagate to uncer- 
tainties in the stellar mass estimates, yielding notably lower masses when stellar age is unknown. 
Given the coarse grid of stellar metallicities, the age-metallicity degeneracy is likely to contribute 
to the large discrepancies in spectroscopic and photometric measurements. Finally, we proposed 
a mass correction to our sample that incorporates these systematic offsets as a function of their 
photometrically determined age and metallicity, finding that the sample median mass increases by 
a factor of roughly 1.3. We would like to comment that these results apply to optical and near- 
infrared broadband photometry of low redshift ETGs. Detailed quantification of the systematic 
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errors involved in SED fitting have produced results consistent with what other studies have found 
and have allowed us to obtain consistent stellar mass estimates for this sample. 

We have used data from the 4th data release of the UKID SS survey, which is d escribed in 
detail in lWarren et^ (12003). The UKIDSS proje ct is defined ir i iLawrence et al.1 (|2007l ). UKIDSS 
uses the UKIR T Wide Field Came ra (WFCAM; ICasah et all (J2007l )). T he photometric system 
is described in iHewett et al.l (J2006l ). and the calibration is described in iHodgkin et al.l (|2009l ). 
The pipeline proces sing and science archive are described in Irwin et al. (2011, in prep) and 
Hamblv et al.l (|2008l ). Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan 
Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department 
of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the 
Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site 
is http : //www . sdss . org. The SDSS is managed by the Astrophysical Research Consortium for the 
Participating Institutions. The Participating Institutions are the American Museum of Natural His- 
tory, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western 
Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced 
Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear 
Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist 
Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max- 
Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New 
Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, 
Princeton University, the United States Naval Observatory, and the University of Washington. The 
authors would also like to thank Olivier Ilbert and Stephane Arnouts for help with using LePhare. 



Table 4: Summary of results: optical+NIR. Offset is the median of the A(log Al) histograms and 2a is the zb2-sigma uncertainty 
on the median. Offset is provided for the fixed— free parameter tests (Figure [TT]) at the median spectroscopic values of t, Z, 
and Ay for BC03 and CBIO models. Unless otherwise noted, these tests assume a spectroscopic redshift, Chabrier IMF, and 
Cardelli extinction law. Note that galaxies fit with the median spectroscopic age of 4.0 Gyr (CBIO) suffered catastrophic fitting 
failures, and so the offset and uncertainties are reported here for an age of 4.5 Gyr. 





Offset 


BC03 

2a 


Offset 


CBIO 

2a 


Note 




spectro-z — photo-z 

BC03 - CBIO 

BC03 (CBIO) - PEGASE.2 


0.0128 

-4.64E-3 

0.0962 


+0.804/-0.500 
+0.118/-0.163 
+0.0837/-0.112 


0.0555 

-4.64E-3 

0.0858 


+0.771/-0.524 
+0.118/-0.163 
+0.0934/-0.132 


0.000 <Zp< 0.300, 9 < log M/M^ 
STELIB spectral library 


3 < 12 

IsD 


Chabrier — Salpeter 


-0.227 


+0.0855/-0.0323 


-0.225 


-h0.0441/-0.0281 


STELIB spectral library 




Cardelli — Calzetti 
Cardelh - CF2000 


-4.50E-3 
-0.104 


+0.0711/-0.140 
+0.142/-0.170 


-l.OOE-4 
-0.0680 


+0.109/-0.0375 
+0.156/-0.131 


Cardelh/Calzetti from LePhare < 
TV = 1.0, n = 0.3 (in GALAXEV 


::ode 
code) 


u— no-n 
ugrizYJHK — ugriz 


-0.0257 
-0.0566 


+0.131/-0.151 
+0.0868/-0.125 


-6.30E-3 
-0.0419 


+0.115/-0.121 
+0.147/-0.108 








Fixed— Free 
t 
Z 

Av 



0.0833 -F0.108/-0.144 
0.0587 -h0.191/-0.0442 
0.0276 -h0.280/-0.147 



0.0590 -F0.284/-0.201 

-1.50E-3 -h0.0273/-0.151 

-0.164 -F0.554/-0.488 



tmed = 3.25 (4.50) Gyr, BC03 (CBIO) 

[Fe/B]rned = 0.0932 (0.0932), BC03 (CBIO) 

Av,med = 0.12 (0.10), BC03 (CBIO) 



Table 5: Summary of results: complete. Offset is the median of the A(log A^) histograms and 2a is the ib2-sigma uncertainty 
on the median. Offset is provided for the fixed— free parameter tests (Figure [TT]) at the median spectroscopic values of t, Z, 
and Ay for BC03 and CBIO models. Unless otherwise noted, these tests assume a spectroscopic redshift, Chabrier IMF, and 
Cardelli extinction law. Again, the reader is cautioned on use of the results in this table, per the reasons outlined in Section [71 







Offset 


BC03 

2a 


Offset 


CBIO 

2a 


Note 


spectro-z 
BC03 - 
BC03 (CBIO) 


- photo-z 

- CBIO 

- PEGASE.2 


0.0281 
0.0131 
0.0871 


+0.780/-0.330 
+0.0424/-0.0326 
+0.0516/-0.0908 


0.0185 
0.0131 
0.0802 


+0.772/-0.316 
+0.0424/-0.0326 
+0.0550/-0.151 


0.000 <zp< 0.300, 9 < logM/Me < 12 
STELIB spectral library 


Chabrier - 


- Salpeter 


-0.233 


+0.0144/-0.0140 


-0.231 


+0.0173/-0.0151 


STELIB spectral library 


Cardelli - 
Cardelh - 


- Calzetti 

- CF2000 


-0.0193 
-0.0693 


+0.0777/-0.0386 
+0.0918/-0.0501 


5.10E-3 
-0.0225 


+0.0578/-0.0718 
+0.0941/-0.0433 


Cardelli/Calzetti from LePhare code 
Ty = 1.0, n = 0.3 (in GALAXEV code) 


U— 1 


lO-U 


-0.0128 


+0.0410/-0.102 


-0.0166 


+0.0571/-0.187 


— 



I 

IsD 

i 
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A. SDSS fiber vs. Kron radius 

In this paper, we provide spectroscopic parameters measured within the SDSS fiber radius of 
1.5". However, photometric magnitudes are measured using an adaptive aperture radius of 3 x rK,i 
(see Section [2]). Figure [T71 shows the difference in flux ratio of SDSS FiberMags to Kron magnitudes 
between the wavebands g and x, with x = riz, as a function of Mr for the complete sample. The 
plot shows that there is essentially no color gradient between the 1.5" and 3 x rK,i apertures, from 
g through z. Indeed, the FiberMag/(3 x rK,i) flux ratios have color consistent with zero inside the 
2(7 error bars for most of the absolute magnitude range, with no noticeable dependence on M^. We 
conclude that it is reasonable to compare STARLIGHT spectroscopic parameters, measured within 
the fiber aperture, to those from SED fitting, measured within the 3 x rK,i aperture. 



B. The initial mass function (IMF) 

Expressions for the initial distribution of stars in a stellar population as a function of mass, 
known as the stellar initial mass functio n (IMF), have been proposed by several authors (e.g. 



Salpeterlll955l : IScaldll986l : IChabrierl l2003l ) . In this paper, we primarily use the Chabrier IMF for 



BC03/CB10, parameterized as 



exp[- ^'°g"^-^°g"^-) ] , for 0.1 < m/Me < 1, 
m~^-^ ,for 1 < m/MQ < 100, 



A/ N I exp — ^ TT^ —I , ior u.i ^ in jvip) \ i, /-n-,\ 



with 771c = O.O8A^0 and a = 0.69, and the Scalo IMF for PEGASE.2, approximated here by two 
power-law segments: 



m ^-^^ ,for 0.1 < m/TW© < 1, 
m-2-67 Jor 1< m/MQ < 120. 






The logarithmic slope of ^(w,) ex m~'^'^^ (B3) for the Salpeter IMF is tested against the Chabrier 
form in Section 17.31 We can compare the results to the analytical approximations, derived by 
integrating Equations (Bl), (B2), and (B3) as m x S,{m). Salpeter masses are expected to be a 
factor of 2.15 (~0.33 dex) larger than Chabrier, and Scalo masses are expected to be larger by a 
factor of 1.05 over Chabrier. The analytical result (0.33 dex) falls within the 3a confidence limits of 
the Chabrier/Salpeter comparison in Figures [5] and [UJ We then take the difference in stellar mass 
between the Chabrier and Scalo IMFs computed above and apply this correction to A(logA4) for 
BC03(CB10)-PEGASE.2, as described in the text. 
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Fig. 1. — Distributions of models, ages, metallicities, and color excesses, using BC03 (left), CBIO 
(center), and PEGASE.2 (right). See Table [6] for "Stellar Population Model ^" information. 
Gray histograms use the spectro-z, and black histograms use the photo-z. Note that PEGASE.2 
templates use solar metallicity. Stellar age is measured using the median of exp(— x^/2) from all 
fits, whereas LePhare only provides model (metallicity) and E(B — V) from the best x^ value. 
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Table 6: Stellar population model parameters (BC03/CB10/PEGASE.2). 





BC03/CB10 


PEGASE.2 


Number 


r 
(Gyr) 


Z 


SFR 


1 


0.1 


0.2 X Zq 


b(t = 0) 


2 


0.1 


0.4 X Zq 


r = 0.1 


3 


0.1 


Zo 


r = 0.3 


4 


0.1 


3x Zq 


r = 0.5 


5 


0.3 


0.2 X Zq 


r = 0.7 


6 


0.3 


0.4 X Zq 


r = 1.0 


7 


0.3 


Ze 


r = 2.0 


8 


0.3 


3x Zq 


r = 3.0 


9 


1.0 


0.2 X Zq 


r = 5.0 


10 


1.0 


0.4 X Zq 


r = 7.0 


11 


1.0 


Zq 


r = 9.0 


12 


1.0 


3x Z© 


T = 10.0 


13 


2.0 


0.2 X Zq 


T = 15.0 


14 


2.0 


0.4 X Zq 


T = 20.0 


15 


2.0 


Zq 


0.5 X lO-^XoMyr^i 
for t < 20 Gyr 


16 


2.0 


'ixZQ 


— 


17 


3.0 


0.2 X Zq 


— 


18 


3.0 


0.4 X Zq 


— 


19 


3.0 


Zq 


— 


20 


3.0 


3x Zq 


— 


21 


5.0 


0.2 X Z© 


— 


22 
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0.4 X Z© 


— 


23 
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Zo 


— 


24 
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3x Z© 


— 


25 


10.0 


0.2 X Zq 


— 


26 
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0.4 X Z© 


— 


27 
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Zq 
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28 
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3 X Z0 
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29 
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0.2 X Zq 
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Fig. 2. — Comparison between the spectroscopic and photometric redshifts, having corrected for 
zero-point magnitude offsets (see Section I4.1.ip with a bright subsample 14.0 < mi < 16.4. Top 
panels show redshifts for the optical+NIR sample, with squares representing the median values in 
each bin (200 objects) and their 2a uncertainties. Bottom panels show residuals, Zg — Zp, normalized 
by the spectro-z. 
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Fig. 3. — Comparison between spectroscopic (gray) and photometric (black) ages, metallicities, 
and extinctions for the optical-|-NIR sample. Stellar population models assume single bursts at 
t = and the parameters in Tabled) STARLIGHT measured extinction has been binned to the six 
E{B — V) values in the range — 0.5 (converted to extinction in magnitudes using Ry = 3.1. Solid 
blue (green) lines show the sample median, measured with STARLIGHT (LePhare), with their 2a 
uncertainties shown as dashed lines. 
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Fig. 4. — Each plot shows differences (STARLIGHT — LePhare) in spectroscopic and photometric 
parameters t, Z, and Ay, as a function of their spectroscopic values, measured with STARLIGHT. 
Black symbols represent median values in each bin. Linear fits {dashed lines) to the data are 
presented, with their 1-sigma uncertainties in the gray shaded region. 
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Fig. 5. — Signal-to-noise ratios (S/N) for the optical+NIR sample. The S/N in the li-band is too 
low to allow the fitting of galaxy images with Sersic models (as in previous papers in this series; see 
Section r2.2p . so SExtractor was used to obtain the magnitudes and corresponding errors. S/N in all 
of the other bands is remarkably high. Dashed lines represent the median S/N in each bandpass. 
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Fig. 6. — Errors in the stellar mass computed from the observed apparent magnitudes minus the 
Monte-Carlo magnitudes (see Section [7. ip . using the spectro-z (gray) and photo-z {black), for the 
optical+NIR sample. The ugrizYJHK magnitudes were varied within their la errors according to 
a normal distribution, and the resulting magnitudes were used to scale the template SED to obtain 
the stellar mass estimate. The solid and dashed lines indicate median and ±2(7 limits, respectively, 
for the gray histogram. 
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Fig. 7. — Histograms of A{logAi)s-p, defined in Section [U for spectroscopic and photometric 
redshifts, using the BC03, CBIO, and PEGASE.2 stellar population models. The solid and dashed 
lines indicate median and ±2(7 limits, respectively. Only objects in the range 9 < log A^/A40 < 12 
are included in this figure. 
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Fig. 8. — Histograms of A(logA^), defined in Section [U for tlie optical+NIR sample and various 
BC03 model parameters. With the exception of the spectro-z— photo-z mass plot {gray), the redshift 
is fixed in each case. A Chabrier IMF and Cardelli extinction law is assumed, unless otherwise 
noted. Overplotted are a vertical dashed line and a Gaussian with a = 0.1 dex centered on zero. 
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Fig. 9. — Histograms of A(logA^), defined in Section [6l for the optical+NIR sample and various 
CBIO model parameters. With the exception of the spectro-z— photo-z mass plot {gray), the redshift 
is fixed in each case. A Chabrier IMF and Cardelli extinction law is assumed, unless otherwise 
noted. Overplotted are a vertical dashed line and a Gaussian with a = 0.1 dex centered on zero. 
A color version of this figure is available online. 
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Fig. 10. — Effect of using different optical+NIR bandpasses for scaling tlie SED (using BC03). 
With known redshift, tliis histogram suggests that there is no systematic bias in choosing all bands 
as opposed to X-band to scale the SED. The median difference (solid vertical) is indistinguishable 
from zero at the O.So" level (dashed vertical). Note that LePhare does not accept only 1-band for 
scaling when the redshift is unknown. 
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Fig. 11. — Plots showing log ^A fixed — logA^/ree offsets, described in Section [U as a function of 
spectroscopically determined age, metallicity, and extinction for the optical-l-NIR sample, using 
M09 (left), BC03 {center), and CBIO (right) models. The color code is ~ from red to blue - # 
galaxies < 50, < 100, < 200, < 300, < 400, < 500, < 600 (black > 600). Measurements in a given 
bin are only considered if the bin contains >40 galaxies without catastrophic fitting failures (as 
returned by LePhare). Overplotted are linear fits to each data set. The error bars denote the 
95% CL on the median, normalized by the square root of the number of galaxies in the given bin. 
Provided within each plot are the mean offset (//), the standard deviation of the residuals (cTres)) 
and the slope and y-intercept for a model line fit from a x^-minimization. 
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Fig. 12. — The mass correction, e, plots as a function of photometrically determined age and 
metallicity, using M09 (left), BC03 (center), and CBIO (right) models. The color code is defined 
in Figure [TTJ Measurements in a given bin are only considered if the bin contains >20 galaxies 
without catastrophic fitting failures (as returned by LePhare). Overplotted are linear fits to each 
data set. The error bars denote the 95% CL on the median, normalized by the square root of 
the number of galaxies in the given bin. Provided within each plot are the mean offset (fi), the 
standard deviation of the residuals (cres), and the slope and y-intercept for a model line fit from a 
X^-minimization. 
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Fig. 13. — Stellar mass distributions for the optical+NIR sample. Only the spectroscopic redshift 
is constrained in the gray histogram. The black shaded (outlined) histograms show the stellar 
masses, corrected for age (metallicity) according to the procedure described in Section [7.8.21 These 
corrections utilize the spectroscopic measurements from M09 {left), BC03 (center), and CBIO 
(right) stellar population models, respectively. 
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Fig. 14. — Best-fit output spectra from LePhare for a single galaxy with age fixed to the spectro- 
scopic value of 12.5 Gyr (solid blue) and photometrically fit to 2.2 Gyr (dash-dot red), selected from 
the optical+NIR sample. The observed apparent magnitudes, converted in flux, are overplotted 
(^squares) with their la errors. The ugrizYJHK filter curves are overplotted for reference. 
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Fig. 15. — Comparison of total stellar masses for the complete sample and the 

same sample of galaxies from a group at the Max Planck Institute (obtained from 
http : //www . mpa-garching . mpg . de/SDSS/DRT/Data/stellarmass . html) , who use SDSS Fiber- 



Mag photometry. The hlach circles are obtained by median-binning the data, the dashed line is 
the 1:1 line, and the dotted line is arbitrarily drawn 0.1 dex below the 1:1 line to separate the two 
trends observed here. 
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Fig. 16. — Distributions of age, star-formation decay time scale, metallicity, and extinction for 
galaxies lying above {top panels) and below (bottom panels) the dotted line in Figure [T5l 
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Fig. 17. — Difference in flux ratio of SDSS FiberMags over Kron magnitudes between the g and x 
wavebands, with x = [riz], as a function of the absolute r-band magnitude for the complete sample 
of ETGs. Black circles are obtained by median-binning the data, with error bars marking the 2cr 
uncertainty on median values. Notice that there is no variation in the difference of flux ratios from 
g through z, i.e. no significant color gradient between the Fiber and Kron apertures. In ^r — z there 
is a small negative offset, but it is less than 0.005 magnitudes. 
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